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Outline of Presentation 



• Silicates are truly the ‘ubiquitous material’ — found everywhere! 

- Natural systems: found in many minerals and rocks 

- Technology — coatings, structural ceramics 


• Apply Knudsen Effusion Mass Spectrometry (KEMS) to study thermochemistry of 
silicates 

- KEMS allows measurement of equilibrium vapor pressures above condensed 
phase 

- Unique challenges 

• Complex vaporization behavior 

• Kinetic barriers to vaporization 


• Examples 

- Geology: Thermochemistry of Olivine 

- Technology: Thermochemistry Y 2 0 3 -Si0 2 and Yb 2 0 3 -Si0 2 coating systems 

• Vaporization kinetics of silicates 
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Silicates in Geology, Mineralogy and Planetary Science 



- Over 90% of the Earth’s crust consists of silicate minerals 

- Olivine (Fe x Mg 1 . x ) 2 Si0 4 primary constituent of earth’s mantle 



39% 



Plagioclase 

Alkali feldspar 

Quartz 

Pyroxene 

Amphibole 

Mica 

Clay 

Other Silicates 
Nonsilicates 


D. Perkins, Mineralogy, 3, Prentice Hall, 2011. 



- Moon, Mars, Asteroids, Comets, Interplanetary dust particles and... 

..Hot, rocky 
exoplanets 
(maybe!) 



G. E. Brown, Rev. Mineral. Geochem., 5, 275-381, 1980. M. T. DeAngelis etal., Am. Mineral., 97, 653-656, 2012. 
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Silicates in Materials Science 

- High-Temperature Materials: Silicate Coatings 




Intake 


Compression 


Combustion 


Exhaust 


Cold Section 


Hot Section 


- Silicon-based ceramics: combustion chambers, static parts in hot stage 

- Protective coating against water vapor, condensed phase deposits 

- Rare earth (RE) silicates (RE 2 0 3 ).n(Si0 2 ) 
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Knudsen Effusion Mass Spectrometry (KEMS) 



Knudsen Cell: 
Condensed Phase/Gas 
Equilibrium 



'' v — Parameter 


- Pyrometer 
sighting pod 


Direct Molecular Beam from 
Effusate 

Into Mass Spectrometer 



Knudsen Cell: 1909 

Couple to mass Spectrometer: 1950s (Ingrahm et al.) 
Continuing valuable applications of these methods! 


www.nasa.gov 



National Aeronautics and Space Administration 

Standard Calibration Material: Au 


• Temperature Measurement: very critical 

• Use pyrometer (non-contact) 

• At triple point: determine calibration constant 


p - 

r M 


SIT 


G 
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Standard Calibration Material: Au 
Heat of Vaporization Checks Temperature Calibration and 

Instrument Response 

A V H° = -R*(-41 .162) = 342.20 kj/mol 



Au (s) = Au (g) 


Tables = 342 kJ/mol 


A V G = A V H -TA V S = -RT In K v = -RT ln( P M ) 
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Mass Spectrometer P M = 


P M = partial pressure of M; 


kIT 


(7 


k = instrument constant; I = ion intenisty; 
T = Absolute temperature; 
a = ionization cross section 
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Proceed to Oxide Solutions 

• Partial pressures <-» activities 

• Complex vaporization and ionization behavior 

• Si0 2 (s) = Si0 2 (g) 

• Si0 2 (g) + e — > Si0 2 + + 2e _ 

-> SiO + + O + 2e _ 

• Si0 2 (s) = SiO (g) + Vfe 0 2 (g) 

• SiO(g) + e ' SiO + + 2e _ 

Si + + O + 2e _ 

• Calculation of cross sections for molecules 

• Vaporization may be kinetically limited 

• Container Issues 

• Need inert container or container with known interactions: Mo, W, Pt, Ir 

• Silicates are very reactive! 
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Solutions: Measure Partial Themodynamic Quantities 

Olivine: Fe0 1 . a (Mg0) 1 . (3 (Si0 2 ) 1 . Y 
Same Phase; Variable Stoichiometry 



Pure Compound : 

FeO(s) = Fe(g) + 1/2 0 2 (g) 



Solution : 

FeO (solution, a < 1) = Fe(g) + 1/2 0 2 (g) 




1/2 


a FeO 



P \P 1 1/2 

r FeL r 0 2 J 

p° r p° il/2 

r Fe L r 0 2 J 


ln(a Fe0 ) vs 1/T — partial molar enthalpy 
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Procedure 

• Ion intensity measurements of relevant species for: 

1 . Pure compound 

2. Solution 

• Best to have in-situ pure compound and solution 




Water cooling 


Precision x-y table — ^ 


x positioner 


y positioner 


CD- *0-63.233 


Knud sen celt — ■ - 
Mo Knudsen-cell holder — - 


W foil -heating element 
Ta heat shields 


Power feed-through 


Pyrometer sight 



• BUT, for the highest temperature (>2000K), need to use one cell and change specimens. 
Assume constant calibration factor. 
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Example I — Natural Systems: Olivine 

• Mg 2 Si04(Forsterite)-Fe 2 Si04 (Fayalite) 

• Primary Constituent of Earth’s Mantle 

• Sources: 

- Pure form found on Hawaii Green Sand Beaches: Volcanic pipeline to Mantle 


• Important in volcanism, meteorites, likely constituent of other planetary bodies 

• Very reactive, particularly above melting. Use Ir cell. 
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Olivine- Starting Material and Characterization 



93% forsterite and 7% fayalite, Fo 93 Fa 7 - (Fe 0 07 Mgo 93 ) 2 Si0 4 



XRD pattern and Rietveld refinement of the as received olivine samples. 


ICP-OES analysis 
of the as received 
olivine samples. 


nent 

*wt (%) 

Al 

0.0120(6) 

Ca 

0.035(2) 

Co 

0.0120(6) 

Cr 

0.052(3) 

Fe 

5.01(3) 

Mg 

30(2) 

Mn 

0.075(4) 

Na 

0.0080(4) 

Ni 

0.27(1) 

Sc 

0.0040(2) 

Si 

20(1) 


^Uncertainties of 
the analyses are 
given in 
parentheses. 


Phase content - 


Forsterite - 87.7 ± 0.3% 
Enstatite - 7.1 ± 0.2% 
Silica - 0.84 ± 0.6% 
Sapphirine - 0.5 ± 0.1% 
.Clinochlore - 3.9 ± 0.2% 


Heating to > 1060°C 
removes impurities 
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Complete van’t Hoff Plot 



T(K) 

2300 2200 2100 2000 1900 1800 1700 



Temperature dependence of ion intensity ratios of 
Mg + , Fe + , SiO + , 0 + and 0 2 + in the olivine sample. 



Fegley and Osborne, Practical Chemical Thermodynamics 
For Geoscientists . Elsevier 2013, Fig. 12-11. 


Measurements show good agreement with the 
phase diagram calculated by Bowen and Shairer. 

Bowen and Schairer, Am. J. Sci. 29, 151-171 (1935). 
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Olivine — Solution of 
Mg 2 Si0 4 (Fosterite)-Fe 2 Si0 4 (Fayalite) 
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Previous Data and Models of Olivine 



• Thermodynamic measurements 

- Nafziger & Muan (1967); Kitayama & Katsura (1968) — from P(0 2 ) 
and stable phases 

- Sakawa et al. (1976): Equilibration method for a(FeO) 

- Plante et al. (1992): KEMS measurements of a(FeO) 

- Wood & Kleppa, Kojitani & Akaogi: Calorimetry 

- General agreement: a(FeO): Positive deviation from ideality 


• Saxena et al. (1993): (Mg,Fe) 2 Si0 4 Regular Solution Lo = 9000 

• Decterov et al.: Sublattice 


• Fabrichnaya (1 998): (Mg,Fe) 2 Si0 4 p and y spinel; subregular solution 
with temperature dependent mixing parameters 
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Thermodynamic Activities in Olivine - (Fe 2 Si 04 ) 0 o 7 (Mg 2 Si 04 )o 9 3 + 

MgO 



Olivine + MgO 
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5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 

1/T(K) * 10 4 


Component 

X 

a (1800K) 

Com pare to 
Ideality 

A/^(1700-1950) 


MgO 

0.62 

0.353 

(-)ve deviations 

30.2 kJ/mol 


“FeO” 

0.047 

0.081 

(+)ve deviations 

212.5 kJ/mol 

Consistent with 
literature 

Si0 2 

0.33 

0.046 

(-)ve deviations 

220.2 kJ/mol 



www.nasa.gov 



National Aeronautics and Space Administration 


Thermodynamic Activities in Olivine - (Fe 2 SiO 4 ) 0 .07(Mg 2 SiO 4 ) 0 93 




















*** ,,< l 

1* I 

Si0 2 

• 

j 

******* 

! 

k 











• * * * 


1 * * # 

* * * i 




[ J 



1 







t 

MgO 










FeO 






















5.05 5.1 5.15 5.2 5.25 5.3 5.35 5.4 5.45 5.5 5.55 

1/T(K) * 10 4 


Component 

X 

a (1900K) 

Com pare to 
ideality 

( 1700 - 1950 ) 


MgO 

0.62 

0.165 

(-)ve 

deviation 

-222.0 kJ/mol 


“FeO” 

0.047 

0.053 

(+)ve 

deviation 

-55.2 kJ/mol 

Consistent with 
literature 

Si0 2 

0.333 

0.341 

(+)ve 

deviation 

116.2 kJ/mol 
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Comparison to Models (FactSage) 



Olivine + 0.01 MgO 
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Thermodynamics of Olivine: Lots to Do! 

• Understand changes in activities on melting 



• Compare partial molar enthalpies to total excess free energy 


• Refine current models: our data suggests some components far 
from ideality 
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Example 2: Rare Earth Silicates 
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SiC and SiC based Composites for 

Heat Engines 

• Strength retained to higher temperatures than metals 

• Lighter weight 

• Fiber Reinforced composites give some fracture toughness 

• Protected by Si0 2 scale 

- Slow growing, good in pure oxygen 

- BUT... Attacked by basic molten salts; volatilized by water 



SiC/SiC CMC HPBR Paralinear 

(1100 °-1300°C, 6 atm; Robinson/Smialek 1998) 
Si(OH) 4 volatility (Opila et al., 1998-2006) 
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Combine Desirable Mechanical Properties of SiC with 
Chemical Inertness of Refractory Oxide 



• Lower activity of silica =^> less reaction 

• Molten salt reaction 

- Na 2 0(s) + Si0 2 (s) = Na 2 0-xSi0 2 



Yb 2 Si0 5 



• Water vapor enhanced volatilization 

- SiC + 3/2 0 2 (g) = Si0 2 + CO(g) 
SiO, + 2 H 2 0(g) = Si(OH) 4 (g) 

- P[Si(OH) 4 ] = K a(Si0 2 ) [P(H 2 0)] 2 


SiC/SiC 


100 pm 


H 2 0(g) ^ ^ S i(OH) 4 (g) I, MOH(g) i 

SiO o. MO 

(Underline indicates in solution) 


Meschter et al., Annu Rev Mater Res 43, 559 (2013) 
N. S. Jacobson, J Am Ceram Soc 97, 1959 (2014) 
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Rare Earth Silicates: Good CTE Match to SiC 
Calculated Y 2 0 3 -Si0 2 Phase Diagram 



THERMO-CALC (2010.08.10:09.24) : 

DATABASE:USER 

AC(0)=1 , N=1 , P=1 .01325E5; 



THERMOCALC (2010.08.10:11.43) : 

DATABASE: USER 

AC(0)=1, T=1600, P=1.01325E5, N=1; 



MOLE FRACTION SI02 


1 :X(SI02), ACR(SI0 
2:X(SI02 ), ACR(Y2C 


Indirect evidence suggests that the Si0 2 thermodynamic activity is 
lower in the Y 2 03 -Y 2 Si 05 and Y 2 Si0 5 -Y 2 Si 2 0 7 regions 
But there are no direct measurements! 


Fabrichnaya-Seifert Database 
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Issues with Measuring a(Si0 2 ) in RE Silicates 

• Vapor pressure of Si0 2 too low to measure in temperature range of interest 

• Need measurable signal for Si0 2 — use reducing agent to make excess SiO(g). Tried 
several, selected Mo or Ta 

- For a(Si0 2 ) > -0.02 

• Mo(s) + 3 SiOo (soln) = 3SiO(g) + MoO s (g) 

- For a(Si0 2 ) < -0.02 

• 2Ta(s) + 2Si0 2 (soln) = 2SiO(g) + TaO(g) + Ta0 2 (g) 

- Note reducing agent must not change solid phase composition 
• Monosilicates + disilicates +Ta - leads to tantalates 

• Need to account for non-equilibrium vaporization 

• SiO overlaps with C0 2 (m/e = 44) 

- Use LN 2 cold finger for improved pumping 

- Shutter to distinguish vapor from cell and background 

- High resolution instrument (in our dreams...) 

- Gettering pump for C0 2 
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Monosilicate + Disilicate 



Y 2 0 3 -Si0 2 + Y 2 0 3 -2Si0 2 Yb 2 0 3 -Si0 2 + Yb 2 0 3 -2Si0 2 


THERMO-CALC (2010.08.10:09.24) : 
DATABASE :USER 
AC(0)=1 , N=1, P=1.01325E5; 
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Work in Two Phase Region: 
Monosilicate + Disilicate 



Three cells: 

• Au (reference) 

• 3Mo + Y2O3 * 2Si0 2 + Y 2 0 3 * Si0 2 

• 3Mo + Si0 2 

• Mo as powder and cell material 
Mo(s) + 3Si0 2 (soln) = 3SiO(g) + Mo0 3 (g) 

- Compare cells 1 and 2 
Note that cell is part of the thermodynamic system: Best way to overcome container issue! 



K 


[P(SiO)J P(Mo0 3 ) 

[a(Si0 2 )] 3 


Mo + 3 Si0 2 = 3SiO + Mo0 3 


a(Si0 2 ) : 


Cell 2: pure Si0 2 

1 _ f [p 0 (5/Q)fp 0 (MoQ3) l 0 ' 33 


K 


Mo + 3 Si0 2 = 3 SiO + Mo0 3 
a ( m> - { [PfSiOfiPWoO,) 


Cell 3: SiO? in silicate 


■ 0.33 


a(Si0 2 ) 


| [l(SiO)] KMoO,) ]° 33 
\[r(SiO)]l°(Mo0 3 )\ 
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Y 2 0 3 .(Si0 2 ) + Y 2 0 3 .2(Si0 2 ) Yb 2 0 3 .(Si0 2 ) + Yb 2 0 3 .2(Si0 2 ) 

T (K) 

1 750 1 700 1 650 K 1 880 1 860 1 840 1 820 1 800 1 780 




Two Phase Mixture 

a(Si0 2 ), 1650K 

Y 2 0 3 .(Si0 2 ) + Y 2 0 3 .2(Si0 2 ) 

0.281 

Y 2 0 3 .(Si0 2 ) + Y 2 0 3 .2(Si0 2 ) 

0.194 
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XRD after KEMS Measurements of RE Monosilicates + Disilicates + Mo: 
Yttrium monosilicate + disilicate +Mo Ytterbium monosilicate + disilicate + Mo 




Position [°2Theta] (Copper (Cu)) 


“ >u * „ 1 . I„. , 1, 1,1 

nr 
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ll H ill 

, 1 III i III III ,111 .1: 1 ,. 1 

,1 lillll ill lilll III , lilli 

ll, II,, lillll 1,1 II ll, ,, 


Y2 ( Si 04 ) 0; Monoclinic; 04-007-4730 1 

1 1 . , ,11,1. 


i h 1 

1, ||. 

I III . 

1 mill ill . 1 ll ill 1. I.ll 

,1,1 1 1 lillll, 1 1 ill 1 II 1. 


Mo; Cubic; 04-004-8483 

1 1 

Y2 Si2 07; MonodiniCj 00-038-0^0 


II 

II 1 

Ill 1 ll 1 „ 

1 1 



Ln2 Si2 07; Monoclinic; 00-021-1014 

1 1, i 1 

, 1 

II 


I 1 ll 1 M 1 

III ll 



1 , 1 1 

Y2 Si2 07; Monoclinia 00-0|42-(|)167 , . 

1 T 1 1 i 

HE 

ILL, 


1 1 ili ll 

l 1 ill 1 III ill 1 

HI ill 'I' 1 Ml 



Phase 


Y 2 0 3 -(Si0 2 ) 


Y 2 0 3 .2(Si0 2 ) 

Mo 


[634] G34 



Phase wt (%) 

Yb 2 0 3 .(Si0 2 ) 56 

Yb 2 0 3 .2(Si0 2 ) 36 

Mo 8 
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Monosilicate + RE 2 0 3 

Y 2 0 3 -Si0 2 Yb 2 0 3 -Si0 2 



THERMO-CALC (2010.08.10:09.24) : 

DATABASE:USER 

AC(0)=1 , N=1, P=1.01325E5; 
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Monosilicate + RE 2 0 3 


Two cells: 

• Au 

• 3Ta + Y2O3 + Y2O3 * Si0 2 

• Ta as powder and cell material — cell is part of system 
2Ta(s) + 3 SiOo (soln) = 3SiO(g) + TaO(g) + Ta0 2 (g) 

- Using P eq (SiO) and FactSage (free energy minimization) 

- Correction for non-equilibrium vaporization 
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XRD after KEMS Measurements of RE Monosilicates + RE 2 0 3 + Ta: 





No. Visible Ref. Code 

Compound Name 

Chemical Formula 

PQS!tinn.L°.?Thfi1-al (Conner mi'll 

Score Display Color 

Qu... Crystal Sy... 

Space... 

SemiQuant [%} 

RIR 

1 

i m 

icQo 04-007-9716 

Yttrium Oxide 

Y2 03 

70 BB Red 

IA... Cubic 

la-3 

49 

9.260 


2 m 

ig8d 04-004-3017 

Tantalum 

Ta 

42 Lime 

IA... Cubic 

lm-3m 

4 

29.920 


3 m 

■f®D 04-001-0709 

Tantalum Silicon 

Ta3 Si 

17^B Puc - 

P;... Tetragonal 

P42/n 

4 

4.270 

1 

4 m 

i($D 04-005-8528 

Yttrium Oxide Si.. 

. Y2 ( Si 04 ) 0 

52 ^B Blue 

P;... Monoclinic 

C2/c 

41 

1.550 


1 5 B 

ioDd 04-011-3996 

Tantalum 

Ta 

1 2 | | Aqua 

IA... Tetragonal 

P-421 m 

3 

7.980 


Ytterbium monosilicate + Yb 2 0 3 + Ta 



1 Yttartaum 0**1* $4ic*1* • Yt*O(8J0t>) 

i i 24 5% 

to 70 

bL_. 

io 

in 


50 

io io ~ r So * 


Phase 

wt (%) 

Phase 

wt (%) 

V 2 0 3 .(Si0 2 ) 

41 

Yb 2 0 3 .(Si0 2 ) 

24 

y 2 o 3 

49 

Yb 2 0 3 

66 

Ta 

4 

Ta 

2 

Ta 3 Si 

4 

Ta 2 Si 

2 
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-1.0 

-1.5 

-2.0 

-2.5 
-3.0 
-3.5 
-4.0 
-4.5 

5.6 5.8 6.0 6.2 6.4 6.6 6.8x10-4 

1/T, K-1 

AH (Si 02 , 1600 K) = (5200.26)' R-2.303 = 99.57 kJ/mol 


y 2 o 3 + y 2 o 3 .(Sio 2 ) 


1700 


1600 

1500 K 

1 

1 

1 

1 

1 

1 



O KEMS 1-2012 
a KEMS 4-2013 





• — Fabrichnaya-Seifert Optimization 
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1 1 



- 1 . 0 - 


-1.5- 


- 2 . 0 - 


c n 


O) 
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-3.0- 


-3.5- 


Yb 2 0 3 + Yb 2 0 3 .(Si0 2 ) 


1650 


T (K) 

1600 1550 



1500 


■2.5 ■ a — A — A — A — A a A 


-a— a — a 


y = -141 2.60(1 /T)-1 .67 


6.0 6.1 6.2 6.3 6.4 6.5 6.6 6.7 

T'lO^K' 1 ) 

AH (Si 0 2 , 1600 K) = (141 2.60)- R-2.303 = 27.05 kJ/mol 


RE 2 0 3(s 1600 K) + Si0 2(s 1600 k) RE2Si0 5(S; 1600 K) 

AHj = measured in this work 

RE 2 Si0 5(s 1600 iq —> RE 2 Si0 5(s 2 98 k) 

AH 2 = H 1600 k - H 298 K 

RE 2 0 3(s , 298 K) - * RElOsfs, 1600 K) 

ah 3 

Si0 2(s , 298 K) ~ ^ Si0 2(s 1600 K) 

ah 4 

2 RE(s, 298 K) + 3/2 0 2(g , 29 g K| — » RE 2 0 3(s , 298 K) 

ah 5 

Si(s,298 K) + 0 2fe . 298 K) Si0 2(s 29 g K) 

a 

< 

2 RE( S , 298 K) + Si( s ,298 k) + 5/2 0 2(gj 298 K) — >■ RE 2 Si0 5(S; 298 K) 

AH 7= ^Hf :RE7 sio^298 K 


AH fj RE silicate, 298 K (kJ/mol) 


Y 2 0 3 .(Si0 2 ) 

Yb 2 0 3 .(Si0 2 ) 


i i 

KEMS Calorimetry* 

-2907 ± 16 -2868.54 ± 5.34 
-2744 ± 11 -2774.75 ±16.48 


a(Si0 2 ), 1650 K 
0.000804 
0.00298 


*Liang et al. “Enthalpy of formation of rare-earth silicates Y 2 Si0 5 and Yb 2 Si0 5 and N-containing silicate Y 10 (SiO 4 ) 6 N 2 ”, J. 
Mater. Res. 14 [4], 1181-1185. **J. A. Duff, J. Phys. Chem. A 110, 13245 (2006) 
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Vaporization Coefficients 

• Vapor Flux (mole/unit area-unit time) leaving a free surface into a vacuum: 
Described by Hertz-Knudsen-Langmuir (HKL) equation 

Jtmttt x 

• Measured flux-Modified by a factor a: Vaporization Coefficient 

J {measured) = 

- Metals: Generally unity; Oxides 10 _1 to 10 -5 ! 

• Free surface vaporization = Langmuir vaporization 

• Important parameter — relatively little expt’l or theoretical work since 1970s 

- True vapor flux in a deposition processes 

- High temperature material vaporization limit 

- True vapor flux in a geochemical/cosmochemical processes 


Jl/rMRT 
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What Leads to non-unity Vaporization Coefficients? 




Vaporization of silica 

- Si0 2 (s) = Si0 2 (g) 

- Si0 2 (s) = SiO(g) + Va 0 2 (g) 

- Si0 2 (s) = SiO(g) + 0(g) 

Complex process 

- Break apart Si0 4 -2 

- Adsorbed Si0 2 (a), SiO(a), 0 2 (a), 0(a) 

- Desorption to Si0 2 (g), SiO(g), 0 2 (g), 0(g) 

- Break 0-0, Si-0 bonds; make 0=0 double 
bond 


Expect a kinetic barrier => flux reduced from 
equilibrium 
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Importance of Vaporization Coefficient 

• Calculate vapor pressures above a condensed phase oxide: 

- Modify by vaporization coefficient 

• Thermodynamic measurements 

- Implicitly assume that a(A(g), solution) = a(A(g), pure component) 

- A(g) — particular species 

• Measurements of these until 1970s, then relatively little work 

• Important parameter has major effect on vapor pressures 
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Summary 

• Knudsen Effusion Mass Spectrometry 

• Powerful tool for thermodynamic measurements 

• In use for many years; but still very useful particularly for solutions 

• Procedures are system specific 



• Example: Olivine 

• Challenge to find ‘inert’ cell material. Iridium probably the best 

• Treat as solid solution of ‘FeO’, MgO, Si0 2 

• The melting point of the olivine sample was determined by the ion intensity discontinuity to be 
1805 °C 

• Compare to standards and derive thermodynamic activities in solid phase. Appears to be 
significant partial molar heats, deviations from ideality 

• Example: Rare-earth silicates 

• The reduced Si0 2 activity in Rare-earth silicates should limit their reactivity with water vapor 

• Solid State rare earth oxides — activity of Si0 2 

• Need reducing agent to obtain a measurable signal for SiO(g), which in turn relates to 
activity of Si0 2 . Reducing agent must not change solid phase composition. 

• Method and choice of reducing agent depends on particular silicate 

• Vaporization Kinetics: Described by vaporization coefficient 
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